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In this paper we study a general nonlinear Schrodinger equation with a time dependent harmonic 
potential. Despite the lack of traslational invariance we find a symmetry trasformation which, 
up from any solution, produces infinitely many others which are centered on classical trajectories. 
The results presented here imply that, not only the center of mass of the wave-packet satisfies the 
Ehrenfest theorem and is decoupled from the dynamics of the wave-packet, but also the shape of the 
solution is independent of the behaviour of the center of the wave. Our findings have implications 
on the dynamics of Bose-Einstein condensates in magnetic traps. 
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I. INTRODUCTION 

One of the most fruitful concepts of Physics is that 
of symmetries. From high-energy physics to condensed 
matter, symmetries play a central role on our under- 
standing of the world. 

As what concerns classical field theories, a symme- 
try is a transformation which preserves the form of the 
equations. In this case the symmetry can help us in 
many different ways. First of all, me may build solu- 
tions which have the same symmetry as the equation. 
Let us take the two-dimensional spatially homogeneous 
nonlinear Schrodinger (NLS) equation 



idtip{r, t) 



+ |V'P]^(r,t). 



(1) 



This equation is invariant under spatial rotations, and 
therefore we are able to search solutions with the given 

symmetry, ip — ip (^\/ + y^^ {x + iy)^ ■ 

Second and most important, Noether's theorem en- 
sures us that once we have found a certain symmetry in 
our model, it is possible to construct certain quantities, 
often with physical relevance, which will be conserved 
during the evolution. For instance, the invariance of Eq. 
(|l|) under time translations, spatial translations and ro- 
tations, give us seven conserved quantities, which are the 
energy 



the linear moment of the center of mass 



„ d , , 



-iV) 



and the angular momentum of the wavepacket 
L = {-ir X V). 



(2) 
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Finally, the existence of symmetries and their asso- 
ciated conservation laws helps us in the study of other 



properties. In particular, the invariance of Eq. (|^) under 
Galilean transformations allows us to rewrite our equa- 
tions on an inertial frame of reference where the center of 
mass is still. This means that the dynamics of the cen- 
ter of mass does not affect at all the dynamics of other 
properties of the wavepacket. 

In this paper we study a generalization of Eq. (0) 
which lacks traslational invariance. Nevertheless we will 
show a Galilean-like symmetry, which allows us to con- 
struct, from any solution, a continuum of other ones 
which follow different classical trajectories. We will show 
that this symmetry implies a decoupling of the dynamics 
of the center of mass with respect to all other proper- 
ties of the wavepacket. We will also point out some very 
relevant applications of our findings to the dynamics of 
Bose-Einstein condensates and of several optical systems. 



II. THE MODEL 

In this work we will consider the following family of 
nonlinear Schrodinger (NLS) equations with a general 
nonlinear term Gdipl) 



idtil^ir,!) 



1 



A + V{r, t)+Gi\ip\) 



V(r,i). (5) 



We will restrict our interest to the case of a quadratic 
potential V{r,t), i.e. 



1 



Vir,t)^-{r,A{t)r). 



A, 



(6) 



Eq. (||) with potential (||) is an accurate model of many 
physical phenomena. In particular it describes the dy- 
namics of a Bose-Einstein condensate in the mean field 
approximation the propagation of optical beams in 
graded index fibers |Q and the propagation of solitary 
waves in fiber trasmission lines with in-line phase modu- 
lators [|. 
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The nonlinear term, G, may adopt many different 
forms depending on the particular application of Eq. 
The most classical cases are the so called power nonlin- 
earities Gdi/^p) = ijT/'lP which arise in mean-field mod- 
els with different spatial dimensionalities. In Nonlinear 
Optics we also find many versions of the so called sat- 
urable nonhnearities, e.g. G(|i/'P) = ±|^/;p/(l-|-/?|?/;p) as 
well as their Taylor aproximations for small u, GdV'P) = 
±|'!/'p — a|'(/'|^. But the nonlinearity need not be local, 
and in applications to Bose-Einstein condensation one 
finds nonlocal expansions of the atom-atom interaction 
GdV-P) = jK{r- r')|?/>(r')pdr', where the kernel is ei- 
ther radially symmetric K{\r — r'\) [^Q or adopts more 
complex dependencies in the case of dipole-dipole inter- 
actions 1^. These are only a few examples of the many 
forms the nonlinear term may have. 

The description of the dynamics involved in a NLS 
equation is of great interest for applications. However, 
except for the very specific one dimensional case with 
G — ±\ip\'^,A = 0, in which the equation may be inte- 
grated by means of the Inverse Scattering method, noth- 
ing can be said about the structure of the solutions. 
There are other tools such as the moment method which 
give us information about the evolution of relevant inte- 
gral quantities characterizing the solution ||^. In some 
cases, these methods are connected to the conformal in- 
variance of some classes of nonlinear Schrodinger equa- 
tions 1^ but have several limitations: (i) They cannot be 
used to build explicit solutions of the equations and (ii) 
they work exactly only on specific cases. To derive a pro- 
cedure which is valid for more general nonlinear problems 
as the ones we consider here one must use some nontrivial 
approximations 

In this paper we will be able to exploit the behavior of 
Eq. (H) with harmonic potential (^) under spatial trasla- 
tions to provide explicit information on a whole class of 
time dependent problems as will be shown below. 



III. BUILDING NEW SOLUTIONS OF THE NLS 
BY SPATIAL TRASLATIONS. 

A. General case 

Let us consider a solution tp{r,t) of Eqs. (|^)-(|^) sat- 
isfying ip{r,t = 0) = ^(r). Our main result is that given 
any solution, ?/;(r,t), there exists a continuum of other 
solutions which are of the form 



(7) 



being R(t) and d{r, t) appropriate functions to be deter- 
mined later. 

To check this point we proceed by inserting the ansatz 
V'R(r,t) given by Eq. into Eq. (^). Using the fact 
that tp(jj t) is a solution of Eq. (||), we are able to cancel 
several terms on both sides of the equation. If we impose 
that the new function i/)r be also a solution of Eq. (||), 



we reach a solvability condition which is made up of all 
the remaining terms 



dR 



i [dte - iAe + (V6i)2 + (2r - R, A{t)R)] V- 



(8) 



This is a set of partial differential equation for the un- 
known function 9{r,t). Fortunately, it is possible to con- 
struct solutions by choosing a linear phase, 



fit), 



(9) 



together with a trajectory R(i) determined by equations 
of Newton type 



fR 



A{t)R = 0. 



(10) 



By applying the hydrodynamic interpretation of the NLS 
equation Q , the precise for of Eq. leads to a diver- 
genceless velocity field, v = V9 = R, which is responsible 
for the global displacement of the solution. 

Finally, we need a global contribution to the phase, 
f{t), which is determined uniquely from 



dt 



dR dR , ,„ . / N \ 
-,-)-(R,^WR) 



(11) 



This contribution can be calculated for each trajectory, 



fit) 



dR dR 

'dt' It 



iR,A{t)R) 



dt. 



(12) 



Therefore, what we get from Eqs. (|^), (^), ( [lO| ) and 
( p^ ) is a new solution of Eq. (|) which is displaced from 
the initial one. It is remarkable that these explicit time- 
dependent solutions are obtained by spatial traslations 
in a system which is not spatially homogeneous and the 
dynamics is defined by simple, linear ordinary differential 
equations. This behavior is exclusive of the harmonic 
oscillator type potential given by Eq. (^) but it is not 
restricted to any specific form of the nonlinear term or 
any dimensionality of the system. 



B. Evolution of stationary states 

A relevant type of solutions of Eq. (||) are the so called 
solitary waves or stationary solutions, which are of the 
form 



V'(r,i) = 0^(r)e^'^*. 



(13) 



The existence and number of these solutions depends on 
the properties of the nonlinear term. In this paper we 
will assume that the nonlinear term is such that these 
solutions exist, which is in fact the case for most choices 
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of G of physical interest Then we may build from 

( |l3| ) new solutions of the type 



i[l-it+9(r,t)] 



(14) 



In this case the whole of the wavepacket moves following 
a classical orbit, while preserving the shape! This inter- 
esting prediction can be confirmed both experimentally 
and numerically. 




FIG. 1. Trajectories of a solution which is initially sta- 
tionary, and is suddenly displaced and imparted an initial 
velocity. We plot the closed trajectories of the center of 
mass in the symmetric confinement (solid line, uj^ = i-Oy = 1, 
R(0) = (1,0), R,(0) = (0,1)) and in the asymmetric trap 
(dashed line, tuy = 1.2, uj^ = 1,R(0) = (1,0), R(0) = (0, 1)). 
The trajectories have been obtained integrating Eq. (H) nu- 
merically. 

In Fig. |l] we show the evolution of two of such 
wavepackets, first in the symmetric trap (solid line) and 
in the asymmetric trap (dashed line). Such solutions were 
obtained by solving Eq. ^ using a split-step method on 
a Fourier basis with 128 x 128 modes. As our analysis 
predicts, the shape of the wavefunction is preserved up 
to the numerical precision of the computer. 



C. Addition of rotational terms. 



where Lz is the Hermitian operator which represents the 
projection of the angular momentum along the rotation 
axis and is given by 



L^V = -i(r, JVV). 



(16) 



The antisymmetric matrix, J, is the generator of the ro- 
tations around the z axis 



J 




(17) 
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The proof presented in Sec. Ill A is also valid when the 
matrix A(t) is non diagonal. An specific case of physical 
interest arises in Bose-Einstein condensation when the 
trap which confines the atoms rotates. In that case it is 
customary to study the system on the frame of reference 
which moves with the trap, at angular speed Q{t). On 
these coordinates the NLS equation reads 



1 



A + Vir)+Gi\iP\)+nL, 



(15) 



FIG. 2. Trajectories of the center of a solution placed 
initially at a::(0) = l,?y(0) = 1, with i?(0) = in two 
difi'erent situations: (a.) = = 1, SI = 1/2, (b) 
Ldx = l,LOy= 3/2, n = 3/2. 

By repeating the same calculations one arrives to a 
classical equation for the wavepacket center, R(t), with 
an additional term due to the centrifugal force 



— + nj— + = 0. 

at^ at 



(18) 
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These equations form a linear system whose solutions 
are easy to obtain. Specifically, for a two-dimensional 
oscillator and constant frequencies, w^j'^j/ and f2, the so- 
lutions are quasiperiodic with characteristic frequencies 
given by 



H 



(19) 



(20) 



Typical solutions are plotted in Fig. ^, where the tra- 
jectory of the wavepacket has been integrated numeri- 
cally up from Eq. ([l8|). It is important to stress the sta- 
bility of these solutions: even in the case of overcritical 
rotation {Vl > ojx.y), when the centrifugal force exceeds 
the restoring force of the harmonic potential, the motion 
is made of bounded oscillations. 

Incidentally, there is a formal equivalence between a 
NLS equation with a rotating trap (15) with lu^ = tOy — 
and a Ginzburg-Landau equation 



idti; = i(-zV + A)V + IV'lV, 



(21) 



with uniform magnetic field. A* = fl{y, —x,0). In this 
model it is particularly intuitive that the wavepacket 
should rotate around the origin, due to the action of the 
uniform magnetic field, just as the above discussed sym- 
metry (O) reveals. 



IV. DECOUPLING THE DYNAMICS OF THE 
CENTER OF MASS 

Up to now we have shown that given a solution '0(r, i), 
we can build many others, ipnij, t), by spatial traslations 
of the initial data. The process can be reversed, so that 
given a wavepacket (p = ip-B.{r,t) which is a solution of 
Eq. (^ we can extract the dynamics of the center of 
mass, R(t), and the internal dynamics of the wavepacket, 
V'(r,t). 

The practical process is as follows. Let 0(r, t) be any 
solution of the NLS equation with a harmonic potential 
(|). The center of mass position is defined as 



Rc(0 = (r)^ J r\^{r,t)\'d" 



(22) 



The dynamics of the center of mass, and of its associated 
momentum (^), is given by Ehrenfest's equations. Us- 
ing the notation from Quantum Mechanics, the expected 
value of an operator A evolves according to 



(24) 



Applying Eq. ( p3| ) to r and to (— iV), we obtain the 
following coupled ordinary differential equations 



^Pc = (-VU) = AR, 
at 



(25) 
(26) 



With some manipulations it is easy to rewrite this system 
as a second order differential equation 



with initial conditions 



AR, = 0, 



dRr 



dt 



Rc(0) = J r|^(r)pd"r, 



t=o 



(27) 

(28) 
(29) 



where ^(r) = 0(r,O) is the initial data of Eq. (||). 

This means that the center of mass already satisfies 
the equations for a valid displacement in our symmetry 
transformation (^. Hence we can define a second wave- 
function, iplr^t), which moves with the center of mass, 
and which is the solution of Eq. (|^) with initial data 



V'(r, 0) = ^(r + Rc) exp ( — ir. 



d-Rc 

dt 



(30) 



This second wavefunction, -(/^(r, t), is located on the center 
of mass 



r\ij{r, i)|^d"r = 0, 



(31) 



and it is the one that carries the dynamics of all observ- 
ables — widths, angular momentum, circulation, etc — 
completely free from the influence of the center of mass. 

Summing up, what all these transformations tell us 
is that if we displace the initial data, or impart some 
speed to its center, we obtain the same solution, ^{r,t), 
centered on different trajectories. 

This result has been obtained with the help of the Er- 
henfest theorem, which states that the center of mass 
should satisfy an equation of Newton type, and which 
was already known [l0|] . However, the result summa- 
rized in Eq. (^) is much stronger since it states that the 
wavepacket is not affected by the dynamics of its centrum, 
as this dynamics can be integrated out of the equations. 



V. APPLICATION TO THE DYNAMICS OF THE 
CENTER OF MASS IN BOSE-EINSTEIN 
CONDENSATES. 



dt 



A={z[Hii,),A]), 



where H{ip) is a nonlinear operator given by 



(23) Ever since the first works with dilute Bose-Einstein 

condensates, there has been an amazingly precise agree- 
ment between theory and experiments. From the studies 
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of normal modes, to the nucleation of vortices, it is usual 
to obtain a good quantitative matching between the pre- 
dictions (let it be collective frequencies or critical speeds) 
and the actual measurements. 

This is most intriguing in the case of experiments which 
involve a mechanical perturbation of the condensate. We 
first focus on the study of the collective excitations of a 
condensate. Such experiments consist of a periodic mod- 
ulation of the confinement of the condensate, and the 
subsequent study of the oscillations of the wavepacket's 
widths. These manipulations have been shown to not 
only modulate the widths, but to induce an exact, and 
extremely strong resonance of the center of mass ]lO| , pl| . 
Nevertheless, both in the experimental results and in 
some rough models, the widths and the center of mass 
seem to be decoupled, thus allowing us to precisely char- 
acterize the normal modes of the condensate. That ob- 
served behavior is easy to understand in the framework 
of the dynamics of displaced solutions described here. 

Another important application of Eq. (|l^ is the study 
of the center of mass of the condensate in the regime of 
overcritical rotation, Q > minjwa;, Wy}. In this regime 
the rotating condensate, which is ruled by Eq. jl^), suf- 
fers a centrifugal force which is stronger than the restor- 
ing force due to the harmonic potential. It is clear that 
in this regime the condensate should be, and in fact it is 
found to be untrapped. 

However, the analysis of the eigenvalues of Eqs. (18) 



which are given by ( [l9| ) proves that the equilibrium point 
at a; = y = is a center and thus dynamically stable. 
Therefore, the only source of instability for the conden- 
sate under overcritical rotations can be due to deforma- 
tions of the cloud. 

This result is a bit more general than the one in , 
where it is proposed the existence of some configurations 
for the condensate, which correspond to centered and el- 
liptically deformed clouds that survive to the action of 
the centrifugal motion. These configurations are stable 
under dipolar perturbations (displacements of the cloud) 
and under quadrupolar excitations (certain type of defor- 
mations). It remains an open problem to show whether 
such states exist which are dynamically stable under any 
deformation. 

As a side result which can be verified in experiments, a 
perturbed condensate in a rotating trap suffers bounded 
oscillations around the origin with two different frequen- 
cies, LUi and UJ2- These frequencies bear a nontrivial de- 
pendency with respect to the angular speed of the trap 
(|l9|), which can be used to better calibrate experiments. 
Finally we must remark that the existence of two differ- 
ent oscillation frequencies for the center of mass, uJi and 
W2, even in the symmetric trap {ux = ^y) represents a 
splitting of the dipolar mode, which is intuitively simi- 
lar to the splitting of the quadrupolar mode due to the 
presence of a vortex. 



VI. CONCLUSIONS AND DISCUSSION. 

In this paper we have built new solutions by simple 
time dependent traslations in a system without trasla- 
tional symmetry. It is remarkable, and probably a spe- 
cial feature of the harmonic potential that this proce- 
dure works. Specially striking is the case of traslation 
of stationary solutions whose center moves harmonically 
whithout any distortion on the shape of the solution itself 
(only a simple phase appears). 

In relation with the previous finding, we have shown 
that the dynamics of the center of mass is decoupled from 
the dynamics of all other properties of the wavepacket. 
This result stands on other works |l^,|l^. However, the 
contribution of this paper is different and stronger, since 
we show that motion of the center of mass can never 
influence any other properties of the wavepacket, let it be 
a Bose-Einstein condensate or, in a similar row, a solitary 
wave made of light. For all these systems the evolution 
will be essentially the same, no matter the initial position 
and initial velocity of the atomic cloud or solitary wave. 

Our calculations are valid for any type of nonlinearity 
which is symmetric under translations, and which de- 
pends only on the density, \^\. This includes the cu- 
bic nonlinearity for Bose-Einstein condensates, G = ji/'P, 
and most reasonable nonlocal terms |^. This, and the 
fact that our calculations do not depend on the dimen- 
sionality of the system, extends the validity of this work 
to condensates with dipolar interactions, charged conden- 
sates, light in Kerr media and light in saturable media. 

The decoupling of the motion of the center of mass 
has also practical consequences. The invariancc of the 
wavepacket dynamics up to displacements and impulses 
on the initial data, explains why it is actually possible to 
measure the frequencies of the normal modes of a conden- 
sate, even when the center of mass of the condensate is 
known to be exponentially influenced by the changes on 
the trapping potential |]lO|,^. This invariance also ben- 
efits experiments with rotating condensates, as we have 
shown above, and a simple analysis reveals an unexpected 
splitting of the dipolar mode of a condensate. 

The situation is different when other type of poten- 
tials are considered such as stationary pinning potentials 
or any non-harmonic trapping potential, such as some 
polynomial candidates, V{x) oc x*, which are being con- 
sidered in the context of all-optical condensation in very 
elongated traps. The presence of such potentials breaks 
our calculations, as the dynamics of the center of mass 
couples to that of the widths by means of these exter- 
nal agents. We wonder if this will imply some new and 
puzzling dynamics in future experiments. 
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